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Abstract 

A large-N diagrammatic approach is used to study coupled quantum dots in a parallel geometry. 
We show that the Friedel sum rule (FSR) holds at lowest order in a 1 /N expansion for this system, 
thereby suggesting that the ground state is a Fermi liquid. Using the FSR together with the 
dot occupancy, we compute the dot system's conductance. Our finding that the expansion 
indicates the system is a Fermi liquid is in agreement with both prior results based on the Bethe 
ansatz and slave boson mean field theory. 

PACS numbers: 73.63Kv, 72.10 Fk, 71.27.+a, 02.30Ik 



1 



I. INTRODUCTION 



Strong correlations in impurity problems have been of tremendous theoretical 
interest .-lii"— This interest has been spurred by the ability to engineer semiconducting 
quantum dots which are highly tunable.-"— In such realizations of quantum dots, both the 
gate voltages of the quantum dot as well as the tunneling amplitudes between the dots and 
the leads can be tuned. This gives one the ability to study these systems throughout their 
entire parameter space. While the first reported dot systems involved single dots,->^ more 
complicated dot structures can now be fabricated.-"— 

With this ability to engineer multi-quantum dot systems comes the ability to realize 
more exotic forms of Kondo physics. In double dot systems, for example, one can explore 
the competition between the Kondo effect and the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction.—"— Because of the non-perturbative nature of the Kondo effect at low temper- 
atures, one might expect this competition to be non-trivial. We argue that it is in fact so 
for double dots arranged in parallel. 

The geometry that we wish to consider is sketched in Fig. 1. Here two closely spaced, 
single level dots do not directly interact with one another, either via tunneling or via a 
capacitive coupling. They are however connected to one another via hopping processes 
involving the two leads. We will want to consider a situation where a total of two electrons 
sit on the two dots and where the leads are symmetrically coupled to the two dots. In such 
a scenario only one continuum channel of electrons (the even combination of the two leads) 
couples to the two dots. In such a system the Kondo effect and the RKKY interaction are 
antithetical to one another. The Kondo effect will promote many-body singlet formation 
while the RKKY interaction will tend to encourage direct triplet formation between the two 
electrons on the two dots. 

One can then imagine (at least) two possible scenarios for ground state formation as 
one lowers the temperature (see Fig. 2). In the first scenario the RKKY interaction dom- 
inates and at some intermediate temperature binds the two electrons on the two dots into 
a triplet. As the temperature is furthered lowered the Kondo effect then comes into play 
partially screening the effective spin-1 object on the dots down to an effective spin-1/2 (i.e. 
an underscreened Kondo effect). Here the physics is non- Fermi liquid with low temperature 
quantities such as conductance having a logarithmic dependence in temperature. This sce- 
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FIG. 1: A schematic of closely spaced quantum dots arranged in parallel. The bare energy level 
on the dots are e^^^^ • 

nario has been argued to hold in a number of studies including those based on a numerical 
renormalization group (NRG).— i^i^i^ The particular implementation of the NRG there, 
however, has been criticized.— 

A second scenario can be envisioned characterized by a lack of a clear separation of the 
RKKY effect from Kondo physics (i.e. the persistence of higher temperature RKKY-induced 
triplet formation down to temperatures below a putative Kondo temperature). In this 
scenario the low temperature state is an overall many-body singlet. At low temperatures any 
effective binding by an RKKY interaction into a triplet state comes undone. Rather the two 
electrons on the two dots participate in the singlet by screening one another in conjunction 
with the lead electrons. The physics in this scenario is sharply differentiated from the first 
by being Fermi liquid. In particular the deviations in the conductance with temperature go 
as T^. This scenario has been supported by both Bethe ansatz computations^^ii^ as well as 
slave boson mean field theory treatments.— 

It is the purpose of this paper to provide supporting evidence for the second scenario. 
In it we present a. U = oo, treatment of the problem. The method of expansion 
is well known^i^ and has been used in the context of quantum dots/impurities.—"— This 
treatment is distinct from the use of slave boson mean field theory (SBMFT) to argue the 
validity of the second scenario. While SBMFT is itself a technique, it represents some 
(uncontrolled) summation of diagrams. Here instead we consider a 1/iV expansion in 
a systematic fashion. In particular we argue that this expansion is consistent with the 
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FIG. 2: Two possible scenarios for how the ground state forms in doubly occupied double dots in 
parallel as the temperature is lowered. In Scenario 1 the two electrons first form a spin-1 object via 
a ferromagnetic RKKY interaction which is then screened to spin- 1/2 by the one effective channel 
of the leads. In Scenario 2 the RKKY interaction does not dominate the Kondo effect and the two 
electrons on the dots form a many-body singlet with the lead electrons. 



Priedel sum rule - a hallmark of Fermi liquid physics and so indicative of the second scenario. 
(In SBMFT, the Friedel sum rule is automatically satisfied by virtue of the quadratic form 
of the action.) 

The paper is organized as follows. In Section II we briefly describe the double-dot model 
we are interested in studying and the 1/A^ method used. In Section III we present results 
for the dot Green's functions, the system's partition function, and the dot occupancy. In 
Section IV we show that the Friedel sum rule holds and so the system appears to be a Fermi 
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liquid, the central result of this paper. In Section V we compute the conductance using the 
Priedel sum rule and compare it to slave boson mean field. We then draw conclusions in 
Section IV. Details related to computing diagrams in the expansion are relegated to 
appendices. 



II. MODEL OF DOUBLE DOTS IN PARALLEL 

The Anderson like Hamiltonian that we study is given by 



/ + 00 
dxc\^d^ci^ + ^ Via {c\„ad^a + h.c. 

(2.1) 

(TCt a 

The cia specify electrons with spin a living on the two leads, I = L,R. The daa specify 
electrons found on the two dots a = 1,2. Electrons can hop from the leads to dots with 
tunneling strength Via- We suppose that there is no interdot Coulomb repulsion and that 
tunneling between the two dots is negligible. We further consider this model in the limit that 
the strength of the Coulomb repulsion on the two dots, Ua, is taken to oo thus precluding 
any double occupancy. The constraint of no double-occupancy is fulfilled by first introducing 
the slave boson formalism, one slave boson for each dot: 

daa = biUa. (2.2) 

Here f^a is the pseudofermion (denoted later in Feynman diagrams by dashed lines) which 
annihilates one "occupied state" on dot a and b}^ is a bosonic operator (denoted later by 
wavy lines) which creates an empty state on dot a. We subject these new degrees of freedom 
to the constraint, 

Qa^biba + Yf^J'^^ (2-3) 
cr 

which we enforce by adding two Lagrange multipliers Ai and A2 to the Hamiltonian 

5H = ^\^Q^ + ^\2Q2. (2.4) 
We thus understand the partition function of the system to be computed via 

Z= Ad\^e'Pi^^+^^)Tre-^i^+'^.^'Q^). (2.5) 



In treating quantities in this model, we will perform a systematic diagrammatic expansion 
in where is the degeneracy of each dot. 

We suppose the dot-lead couplings, Via, satisfy the following ratio condition, 

^ = ^ = A. (2.6) 

This condition ensures that only one effective channel of lead electrons couples to the two 
dots. To see this we reexpress the Hamiltonian in terms of even and odd channels, 

where Fq = (^1^ + V'g^ct) /2. In terms of these new degrees of freedom, the Hamiltonian divides 
itself into an even and an odd sector, 

H = "He + T^o'-i 

/oo 
dx cl^d^Cea + ^ \/2I\{cl^^d„a + h.c) 

era GL 

/oo 
dx cl^d^Coa, (2.7) 

and we see that the odd sector of the Hamiltonian is decoupled from the dots. 

III. SUMMARY OF COMPUTATIONS 

In this section we will compute the Green's functions of the dot degrees of freedom, the 
dot partition function, and the dot occupancy to leading order and where necessary to prove 
the FSR the subleading order in 

A. Dot Greens functions 

In terms of the slave bosons and fermions, the dot temperature Green's functions are 
defined as 

^ {biir)U{r)fi^{0)b,{0)e~fo'^-(^^'no ,,,, 



where the a, /5 indicate which dot the electron is sitting on and Hy is the part of the 
Hamiltonian couphng the leads to the dots. As computed in the appendices, the lowest 
order contributions in to these Green functions are 
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(3.3) 



In Eqn. (13.21) one can make the replacement 1 ^ 2 to get the corresponding expressions 
for G22^\iuin) and G^i'^\iijJn)- We see that the interdot correlation function, G12, first 
sees a contribution at 0{1/N). Here ek denotes the energy of the conduction electrons and 
fk = g^fc_)_]^ is the Fermi function. Finally, we note that -^01,2 is the most negative solution 
of w — -ReSi 2(w) = and we correspondingly define T^^ 2 T^^ a = ~ -^01,2- 

For the purposes of proving the Friedel sum rule, we will need not only these contributions 
to the dot Green's functions but the imaginary part of the subleading {0{1/N)) correction 
to Gil and G22 as well. These follows: 
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i^f(^oi)ri 



where S{uj) ,R{uj), and T{u) are defined by 



n9{-co - Ta,)R{uj) - Txd{uj - TaJ [S{u) + T{co)] , (3.4) 
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(3.5) 
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We obtain G22^^^^ by swapping 1 o 2. 

The above expression for ImG^j^'^''^'* undergoes a dramatic simplification when oj is set to 
0. We find 

Im 



G°"«(. = 0)l=-^i§^. (3.6) 
This simphfication will be useful for proving the Friedel sum rule. 

B. Dot occupancy 

The partition function and the dot occupancy obey the simple relation 

ld\ogZ 

The partition function at 0{1) is given by 

Z = e-^■^ole-^■^o^ (3.8) 
and -^01/02 are defined as in Section III A. After some algebra we then find that 

ndi,2 = — , (3.9) 

1 + 1^1,2 

where 

/^l,2 = -7f • (3.10) 

7riAi,2 

This 0(1) computation of the dot occupancy is sufficient for proving the Friedel sum rule 
at leading order in 1/A^. 



IV. PROOF OF THE FRIEDEL SUM RULE FOR DOUBLE QUANTUM DOTS 

The Friedel sum rule (FSR) states for a system with N-fold degeneracy that the scattering 
phase, S{E), of an electron at the Fermi surface -E = is given in terms of the total 
occupancy, rid, of the dots by 

^(0) = (4.1) 

The Friedel sum rule implies that the zero temperature linear response conductance is then 
given by 



where A = Via/Voa, a = 1, 2. 



From Ref. 



371 . the hnear response conductance can be expressed in terms of the dot 



Green's functions as follows: 



G 



-Im[IV{riG^}] 



(4.3) 



h (1 + A2) 

where G^ is the full retarded Green's function matrix in dot-space and is the bare 
hybridization matrix. These are defined as 



G' 



and 



F, 
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V^llV^12 



(4.4) 



V^llV^12 V, 
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(4.5) 



Thus to show that the FSR holds we need to show that 



4A2 



sm 



\ N ) 



/i 1 + A2 



/i(l + A2)2 

holds in terms of a expansion. 

The left hand side of (14. 6 p to 0(-^) is equal to 

4A2 . o /vrn^x 4A2 . o /vr (ra^i + 72^2) 



Im[Tr{FiG^}] 



(4.6) 



:i + A 



2^2 



. 2 (T^"^d\ 

sm 



:i + A 

4A2 



2^2 



sm 



TT 



n 



dl 



.0(1: 



n 



d2 



2n 



0(1) 0(1] 



di 



n 



d2 



(4.7) 



(l + A2)2iV2 

where we only consider the dot occupancies to 0{1). The right hand side of (14. 6 p to the 
same order equals 
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-Im[Tr{Fz.G"^}] 



(1 + A2 
+ V^^V^^G 



Iml VlG 



0(1/N) 



ii'-^ii 



^0(1/N) 



o(^''n)^^,.gO(i/n)| 



12'-^22 



(4. 



where all the Green's functions were evaluated to the given order in the previous section. 

^^^^ and so does not appear in the above. Using Eqs. 13.61 and one 



Note that ImG 



11/22 

can then readily verify that the FSR (14. 6 p holds for double dots to leading order in a (9(1/A^) 
expansion. This then matches the proof of the FSR in single dots in a large-N expansion 



performed in Ref. 
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T=0 Conductance: Symmetric Dot-Lead Couplings 
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FIG. 3: Plot showing the hnear response conductance for two different values of separation between 
the dots' chemical potential. In the left panel we present results arrived at from the application of 
SBMFT done in Ref. 
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and in the right panel the results coming from the large-N expansion. 



V. EVALUATION OF CONDUCTANCE USING FSR AND 1/N DIAGRAMS 

Having proven the FSR, we will now use it to compute the linear conductance, G, i.e. we 
will use the expression 

^2 . o / vr [rid, + 



G 



■ sm 



(5.1) 



h{l + A2)2 \ N 
where ra^^ and are calculated from diagrams of the partition function Z (see Section III 
B). 

In Fig. 3 we present the conductance for two different spacings between the chemical 
potentials of the two dots (one large and one small). We compare the large-N computation 



done in this paper with the results of Ref. 



31 



We find remarkable agreement between 



the two approaches. Not only are the general trends in the conductances similar, but more 
detailed features such as inflection points are reproduced. As the dot chemical potentials 
are lowered we see that the conductance vanishes. 

In Fig. 4 we plot the corresponding total dot occupancies both from the 1/A^ expansion 
(to lowest order in 1/A^) and SBMFT (from Ref. l3ll ). We again see that there is a close 
correspondence between the two methods. This is notable as in computing the dot occu- 
pancies from the 1/A^ expansion we only keep 0{1) diagrams while the SBMFT accounts 
for diagrams at all orders of 1/A^. 
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T=0 Displaced Electrons: Symmetric Dot-Lead Couplings 
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FIG. 4: Plot showing the total occupation of the dots for two different values of separation between 
the dots' chemical potential. As with Fig. 2 we show the values of the dot occupancy both computed 
in SBMFT and in the 1/N expansion. 

VI. DISCUSSIONS AND CONCLUSIONS 



In this paper we have studied a double- impurity model using a expansion where 
is the degeneracy of each dot level. Using this 1/A^ expansion we demonstrated that the 
Friedel sum rule holds, so providing complementary evidence that double dots in parallel in 
the absence of interdot and capacitive coupling have Fermi liquid ground states. In previous 
works we have argued that this ground state is Fermi liquid using both the Bethe ansatzi^*^ 
and SBMFT.— This then provides a complementary piece of evidence for this Fermi liquid 
nature. 

Such evidence is needed because there exists in the literature an o ppo site set of results 
coming from numerical renormalization group analyses.— i^i^i^ In Ref. |3l|we have suggested 
a possible reason for this discrepancy, namely in discarding modes that arise from the use 
of logarithmic basis, the NRG discards additional effective screening channels that would 
completely screen any putative spin-1 object that might form under a ferromagnetic RKKY 
interaction. But as also pointed out in Ref. |3l|, another possible reason for the discrepancy 



exists. In treating the double dots in parallel all of these methodologies (Bethe ansatz, 
SBMFT, 1/A^, NRG) treat the dots as point scatterers in a continuum media. Then with 
two dots in close proximity to one another there exists the possibility that the physics 
might be influenced by the particular implementation of a UV cutoff that the methodology 
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employs. If the physics was sensitive to this implementation, one would expect to find non- 
universal results. It is notable then that we were able to conclude in this study that the 
expansion produces results (i.e. a Fermi liquid ground state) equivalent to that of SBMFT 
and the Bethe ansatz. This finding then gives additional weight to the universal nature of 
the Fermi liquid ground state of double dots in parallel. 

We will end on a tangential note by stressing the remarkable similarity between results 
obtained for a double dot in a (low order) expansion and SBMFT. This would not 
be necessarily expected as the SBMFT represents some uncontrolled sum over diagrams at 
all orders in and is merely 2. It is an interesting question why these higher order 
diagrams are suppressed strongly even when N = 2. We note though that while the 
expansion agrees quantitatively with SBMFT and while both of these methods agree on a 
qualitative level with the Bethe ansatz, these methods do not produce fine details of the 
conductance of the double dots brought out by the Bethe ansatz (see Ref. l3l|), presumably 
the most trustworthy method of the three as it represents an exact solution of the problem. 
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12 



Appendix A: 0(1) diagrams of Gii{iujn) 



Following the definition of the Green's function in Eq. 13. ![ Gu, can be written as 

a product of diagrams. These diagrams are given in Figs. M and IH] and these are the only 
diagrams at 0(1) that survive the projection 

cttT 



involved in computing any correlation function. The argument j = 1,2 of Zi{j) and 
Li{j, iujn) indicates the diagram involves dot electrons on dot j. Li{l, iun) involves a dressed 
boson propagator, B{z), for dot 1 with self energy Si - defined in Eq. 13.31 This dressing 
occurs already at 0{1) and is straightforwardly obtained via a Dyson equation (see Fig. 4 



in Ref. 



32 as well as Ref. 



33|): 



where Bo{z) = l/{z — iX) is the bare propagator of the boson line. Li{l, iun) is then equal 
to 

Li(Mc..) = r ^e^^^^-. —4- — ^ ' .(A2) 

P ^ J-T,T 27r lUrn " Ml - Si(zi^m - zAi) l[uJn + I'm) " «Ai - 6^1 

In Eqn. IA2l the summation over Matsubara boson frequencies can be performed by the usual 
technique of converting the summation into a contour integration. Doing so we find 

L^{l.ii^n) = r ^e^'''KEo^+^X,)^^^^. (A3) 

Here Eqi is the most negative solution of 2; — ReT,i{z) = 0, b{z) = (e^^ — 1)~^ is the Bose 
function, and Ki is given by 

ay 

K,{E,,) = [1 - l^E,, . (A4) 

In order to perform the A integration we note that the Bose function can be expanded as 
h{z + zAi,2) = e"^(^+*^i'2) + 0(e-2*^^i.2). From this it follows straightforwardly that 

e-^^o^Ki(Eoi) 

Ll(l,2W„) = — : . (A5) 

lUJn - -LAi 
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ICO, 



Li(14coJ 




L2(l,i0)„) 




L4(ia)J 



1 ia)„ 



FIG. 5: Given above are the diagrams involved in the evaluation of the numerator of the expression 
for Gil at both 0(1) (Li(l,iu;„)) and 0(1 /A^) (L2(l, L3(l, icj^), and L4(ia;,i)). The argument 
j of Li^2,3(j> i^n) indicates the diagram involves only dot electrons of species j. (And so L4 does not 
have this argument as it involves both.) The wavy thick lines indicate dressed boson propagators, 
the dashed lines, bare pseudofermion lines, and the solid (internal) lines, conduction electron 
propagators. The labels and colors ((1, black) and (2, red)) indicate the lines belong to the bosons, 
6, and fermions, /, of different dots. 



To determine £"01, we note that 



Si(^) = 2ivri^ 



fk 



(A6) 



evaluates to 



ReSi(a;) = ^ log 



D 



ImSi(a;) 



{u - edi)iVri. 



(AT) 



Eq. IA7I is obtained after converting the summation into an integration over k which brings 
in a factor of As Eqi is the most negative solution of w — ReSi((X') = which satisfies 
-^'Oi < ^di we can write T^^ = e^^ — Eqi and see that T^j^ is always positive. T/i^ is, in effect, 
the renormalized dot chemical potential of dot 1. 
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We will now compute Zi{l) a bubble diagram (see Fig. [5]). It differs from Li{l) in that it 
also involves a conduction electron line (represented as a solid line in Fig. [5]) and involves no 
external frequency. It is important to note that the boson propagator in this case is different 
from the boson propagator in the diagram represented by Li{l, iun)- The boson propagator 
in Ziil) involves summing over a class of diagrams, each of which comes with an additional 
factor of 1/n.—^ In order to make this sum correspond to the usual geometric series^^i^^ we 
introduce a coupling constant parameter g by the transformation of the coupling constant 
Fi — )■ g^Ti. By suitably differentiating the series by g one can turn this series into a regular 
geometric series.—"— We thus arrive at the expression: 

2Ng'^Ti 1 1 1 , 

w ^ (A8) 

/32 iu^ - iX^ - Si(zz/„ - iXi,g) + w„) - (e^a + «Ai) iUn - ek ' 

After performing the Matsubara summation over a;„ and the projection integral, and trans- 
forming the expression for the Matsubara sum over to a contour integral, we obtain 



1 



dg f dz Si (2,51) 



J g J m z- i:i{z,g) 

o Ta 

Proceeding with the coupling constant integration in Eqn. IA8t we find. 

If we now perform an integration by parts and then compute the contour integration about 
the pole at z = E02, we find, 

Z,{l) = e-^^"K (All) 



We finish up then with 



Gil i^^n) = -. 7^- (A12) 

tUJn - -La. 



A similar expression for G^2^\iuJn) is found by interchanging all species indices 1^2. 



Appendix B: 0{1/N) contribution to Gii{iujn) 

qO(i/n) gjygj^ following set of contributions, 

^o(i/7V) _ (-^2(1, iun) + -^3(1, iUn))Zl{2) + Li(l, iUn)Z2{2) + L4{iuj„ 
~ ^i(l)^i(2) 
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FIG. 6: Given above are the diagrams involved in the evaluation of the denominator of the expres- 
sion for Gil at both 0(1) (Zi(l)) and 0{l/N) (^2(1) and Z3). 



Li(l,2a;„)Zi(2)(Zi(l)Z2(2) + ^1(2)^2(1) + Z3) 
Zi(l)2Zi(2)2 

L2(l,iw„) + L3(l,ia;„) Li{iun) Li(l, za;„)Z2(l) Li(l,za;„)Z3 



(Bl) 



Zi(l) ^i(l)^i(2) Zi(l)2 Zi(l)2Zi(2)- 

The first and third terms form the 0{1/N) contribution to Gu of a single dot while the 
second and fourth terms reflect the presence of a second dot. We consider the latter first. 



1. Evaluation of L4(ia;„) and Z3 

We will first consider the diagrams which contain both dot electron species, Li{iuJn) and 
Z3 (see Figs. 5 and 6). We will show that both of them go to zero in the low temperature 
limit. Let us look at L^^iiUn) first: 

^ ,. , 4iVFiF2 ^ 1 1 



, ,,K^m + OJn) - (Qi + iAi) {iPrn " ?Ai - Si(iz/m - iXl))' 

1 1 1 



X 



^ ^ ]_ (B2) 

K^n + ^'n) - i^d2 + iX2) iv'^ - 1X2 - ^2{i^'ra ' ^^2) ' 

Performing the Matsubara sums over v'^ and w^, we arrive at an expression given by (drop- 
ping terms that will disappear once we apply the projection): 

^ 4iVFiF2 >^ 1 
= -2 2^ - 



X 



fiek) 1 



{ium - iXi - - iXi)y e^. - e^/ iu^ + - {iXi + e^i) 



16 



X 



K{ £'o2)e~^^*'^2+-Eo2) g-/3(ed2-efe+jA2) 



(B3) 



Cfc — Ta2 erf2 — Cfc — ^2\—^k + (-(12) 

Of the two terms in the square brackets in the above expression, the first dominates in the 
low temperature hmit as it is proportional to e~^^°'^ and £'02 < £^2 — Cfc- We can thus neglect 
the second term in this limit. There are other such terms which are subdominant as T goes 
to and which are indicated by the ellipses. 

If we perform the final Matsubara sum together with the k' sum we obtain an expression 
involving terms of the form 

dkf{e,) log -— ) (B4) 



and 



-D 



u 

/ 



It is straightforward to show that these integrals vanish as log(-D)/D in the large bandwidth 
limit. 

Now we will calculate the bubble diagram Z3 (see Fig. 6). We have 

11 1 



Z3 oc 



1 1 1 


iVrn' - i^2 - S2(il^m' - ^^2) «(^'m + ^n) " (Qi + ^Ai) i{ujn + v'^) - (Crfa + iA2) ' 

Performing the Matsubara sums as well as the k' sum leaves us with an expression again 
involving terms of the form found in Eqn. (IBSp and so like Li[iujn)-, goes to zero in the 
infinite bandwidth D limit. 



2. Evaluation of remaining diagrams contributing to G^/^^^^ 
In the zero temperature, large bandwidth limit, G'^^^^^'' reduces to 



^o{i/N) _ L2{l,iu]n) + L-i{l,iUn) Li (1, ia;^) ^2 (1) ,^3^. 



But this is precisely the 0{1/N) correction that is seen in single level dots. We can thus use 



Ref. 



32l to complete the evaluation of 



Oil/N)_ 



G^^^^'^\iUJn) = Q'{iOJn) + Bi{iuJn) + B2{iUJn) 
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->- 



M(icoJ 



FIG. 7: Diagram M{iojn) contributing to Gi2{iuJn) at 0{1/N). 



Q'{iUJn) 



z=Eoi ) 



(^01 + i(^n - Cfc - Si (^01 + i^n - Cfc)) 



01 j 



f{k) 



Ek - TaY {Eqi - iuJn + Cfc - Si(Eoi - iUn + Cfc)) 



+2ivri5^/,(i-/fc 



(efc - TaJ2 (ia;„ + e^. - ey - Ta, 
1 



X ■ 



{Eqi + efc - efc/ - Si(Eoi + efc - tk')Y J' 
We have used the same notation as in Ref. 32 



(B7) 



Appendix C: 0{1/N) diagrams for Gi2{iuJn) 

The lowest order contribution to Gi2{iu!n) occurs at 0{1/N) and can be written as 

= Jp^. (CI) 
^2 ^ ' Zi(l)Zi(2) ^ ' 

We already have Zi(l/2) from Appendix A. The diagram representing M{iun) is given in 
Fig. [71 Using calculations similar to Appendix |X] one easily arrives at 

M(..„) ^ f:!!™i)(-,vr;f;)f:!=!^. (C2) 

lUJn - Ta^ ^^n - Ta^ 

which in combination with the expression for Zi(l/2) gives, 

G°r\«-n) = -!^^(-'^^.)^^^. (C3) 
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Note that the term a/TiT^ in Eq. IC3I results from the conduction hne in Fig. [71 
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